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Abstract 

In the context of the Minimal Supersymmetric extension of the Standard Model, 
we study the production of Higgs bosons at the Large Hadron Collider via cascade 
decays of scalar quarks and gluinos. We focus on the cascades involving heavier 
charginos and neutralinos, which decay into the neutral h, A, H and charged 
bosons and lighter charginos and neutralinos, but we will also discuss direct decays 
of third-generation squarks into their lighter partners and Higgs bosons as well as 
top quark decays into bosons. We show that the production rates of relatively 
light Higgs bosons, M$ < 250 GeV, via these mechanisms can be rather large in 
some areas of the parameter space. Performing a fast detector simulation analysis 
that takes into account the signals and the various backgrounds, we show that the 
detection of the neutral Higgs bosons through their decays into bb pairs, and of the 
charged Higgs particles through the t^v signature, is possible at the LHC. 
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1. Introduction 

The issue of detecting the extended Higgs particle spectrum [1] present in the Minimal 
Supersymmetric Standard Model (MSSM) [2] at hadron colliders has triggered a large 
theoretical and experimental activity in the last twenty years; for recent reviews, see for 
instance Ref. [3]. In particular, thorough experimental analyses and simulations [4] have 
shown that at least the lightest MSSM CP-even Higgs particle h should be found at the 
Large Hadron Collider (LHC). In large areas of the parameter space that have not been 
ruled out by LEP2 searches [5], the heavier neutral CP-even H boson, the CP-odd A 
boson, as well as the charged Higgs particles H ± , can also be discovered for an integrated 
luminosity as high as / £ > 100 fb -1 . 

These conclusions were reached by simply investigating direct Higgs boson production 
through Standard Model (SM)-like processes 1 : mainly the gluon-gluon fusion mechanism 
mediated by heavy-quark loops, gg — > h,H,A [8], and the associated production with 
heavy quarks, gg/qq — > h,H,A + bb or tt [9], for the neutral Higgs particles and top 
quark decays, t — > H + b, or associated production with top quarks gg/qq — > H + bi and 
gb — > H~t [10], for the charged Higgs boson. The production cross sections of most of 
these processes [those involving b quarks] are strongly enhanced for large values of tan /3, 
the ratio of the vacuum expectation values of the two Higgs-doublet fields that break the 
electroweak symmetry in the MSSM. 

Another potential source of MSSM Higgs bosons at the LHC is due to the cascade 
decays of squarks and gluinos: because of their strong interactions, these are copiously 
produced in hadronic collisions. These particles could then decay into the heavy chargino 
and neutralinos xf an d Xs a an d) if enough phase space is available, the latter could decay 
into the lighter chargino and neutralinos Xi an d Xi,2 an d Higgs bosons: 

pp -> gg,qq,qq*,qg -> xf, xl, xl + X 

- xtxlxl + K^A,^ +X (1) 

There is also the possibility of a direct decay of squarks and gluinos into the lightest 
charginos xt an d the next-to-lightest neutralinos xt which decay, again if enough phase 
space is available, into the lightest neutralino and Higgs bosons: 

pp -> gg,qq,qq*,qg -> Xi,X2 + x 
- x? + # ± ,/i,#,A +x (2) 

1 It is also assumed, in general, that the Higgs bosons decay only into SM particles and that loop 
contributions of SUSY particles do not substantially alter the rates of the production and decay processes. 
If these assumptions do not hold, as might be the case if SUSY particles are relatively light, the searches 
could become more complicated; for discussions see for instance Refs. [6, 7]. 
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From now on, we will call the decay chain in eq. (1) the "big cascade" while the one in 
eq. (2) is dubbed the "little cascade". Other possibilities for Higgs particle production 
are the direct decays of heavier top and bottom squarks into the lighter ones and Higgs 
bosons, if large enough squark mass splitting is available [11]: 

pp -> t 2 t*,b 2 b 2 with t 2 (b 2 ) ->fi(6i) + h/H/Aorb 1 (t 1 ) + H ± (3) 

as well as top quarks originating from SUSY particle cascades, decaying into H* 1 bosons: 

pp gg,qq,qq*,qg -> t/t + X -> if* + X (4) 

The production of the lightest h boson from cascade decays of strongly interacting 
SUSY particles into the next-to-lightest neutralinos x 2 , which then decay into the h 
boson and the lightest neutralinos x 2 ~* hXii nas been known for quite some time [12]. 
Cascade decays of squarks and gluinos into relatively light H, A and bosons, where 
the main contribution is due to the cascades xti X 2 ~^ X?+ Higgs bosons, have also been 
discussed in the past; see Ref. [13]. Charged Higgs bosons produced via squark/gluino 
decays have recently been considered [14] in the range where they cannot be produced 
from top decays [i.e. for M#-± > m t ] or in association with top quarks [i.e. in the region 
of moderate tan/5 values where the H + tb Yukawa coupling is not enhanced enough]. In 
the present paper, we extend the previous analyses in three directions: 

- We investigate in detail the cascades of strongly interacting SUSY particles into the 
heavier neutral Higgs particles H and A. In particular, we will analyse the production 
rates in the situation where these particles have relatively small masses, Mh,a ^ 150 GeV, 
and moderate to large Yukawa couplings to 6-quarks [the so-called "intense-coupling 
regime" discussed in Ref. [15] for instance]. We will also investigate in more detail the 
cascades of strongly interacting SUSY particles into the lightest h boson, where both 
types of cascades, as indicated in eqs. (1) and (2), are in general present. 

- In the case of the charged Higgs boson, we extend the analysis of Ref. [14] to 
masses below the top quark mass, where the direct decay process t — > H + b is also at 
work and acts as an important source of "background" events. For the production signal, 
we will not only discuss H ± bosons originating from the decays of heavier gauginos, but 
also the more complicated situation where they also come from decays of top quarks 
originating from the cascade decays of squarks and gluinos as in eq. (4). 

- We perform a full Monte Carlo simulation, in which we consider the signals along 
with the various SM and SUSY backgrounds. Using simple kinematical cuts to reduce the 
most dangerous backgrounds and, to be as realistic as possible, simulating the response of 
the CMS detector [the conclusions are expected to remain valid for ATLAS] , in particular, 
the ability to tag b quarks and r leptons, we show that these processes can indeed be 
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detected in some regions of the MSSM parameter space, particularly in areas where there 
is no access to the heavier Higgs bosons in the direct production mechanisms. 

We will not attempt to make an extremely detailed analysis of this topic. The subject 
being too involved and complicated, it calls for rigorous and dedicated simulation studies 
to be covered completely. Rather, our aim will be to open the Pandora box and investigate 
in various possible scenarios [which are hopefully representative of what could broadly 
happen in the MSSM] this new source of Higgs bosons at the LHC and show that these 
particles can indeed be detected. There are two main reasons, in our opinion, which 
motivate such a study [in addition, of course, to the trivial reason that it is a new source 
of MSSM Higgs bosons of unforeseen potential and as such, it must be analysed anyway]: 

The first reason is that an important issue at the LHC, once supersymmetric particles 
are found, will be the reconstruction of the SUSY Lagrangian at the low-energy scale, 
which would allow the structure of the fundamental theory at high scales to be derived. 
This can be achieved only by measuring some of the couplings between the SUSY and 
the other particles. Among these, the couplings of superparticles to Higgs bosons are of 
special importance, since they also probe the electroweak symmetry breaking sector and 
might decide which Higgs scenario is at work. The cascade decays involve these couplings 
and would provide crucial information to achieve this goal. 

Another argument for this study is that, in the MSSM parameter space, there is a hole 
in the region with 130 GeV < Ma ^ 170 GeV and tan/5 ~ 5, where only the lightest 
h boson can be found at the LHC [a hole also exists for Ma larger than ~ 200 GeV 
and tan/3 ~ 5-10]. This is because the dominant production processes for the heavier 
neutral Higgs bosons, gg — > H/A + bb or tt and gg — > H/A, have not sufficiently large 
cross sections [the Yukawa couplings of the H and A bosons to 6-quarks, proportional to 
tan f3, are not sufficiently enhanced, while the couplings to t-quarks are suppressed] and 
the charged Higgs boson cannot be probed in top decays since M H ± > m t — and its 
coupling to tb pairs [which is also proportional to m h tan (3 and m t / tan 0\ is not sufficiently 
enhanced. We will show that in the cascade decays, the value of tan (3 does not play a 
crucial role and, therefore, this hole could be filled up by searches through these decays. 

In the next section, we discuss the production of squarks and gluinos at the LHC 
and their decays into the MSSM Higgs bosons through cascades involving gauginos. We 
also discuss the direct decays of top and bottom squarks to the Higgs bosons and for 
m t > m H ± , H ± production from direct decays of top quarks coming from q and g decays. 
In all cases, the cross sections times branching ratios will be presented in various scenarios. 
In section 3, we present our event generator analysis in selected areas of the parameter 
space and show that the signal events can be detected in spite of all the SM and SUSY 
background processes. A short conclusion will be presented in section 4. 



4 



2. The signal cross sections 

At the LHC, the total squark and gluino cross section in the various pair and associated 
production processes listed in eqs. (1) and (2) is of the order of a(q + g) ~ 110 to 3 
pb for sparticle masses m g ~ rrig ~ 0.5 to 1 TeV, leading to a large, ~ 3x 10 7 to 10 6 , 
number of events with an accumulated luminosity of J C ~ 300 fb -1 . These squarks and 
gluinos can decay into the heavier chargino and neutralinos X21X3 an d xli with large 
branching fractions of about a few times ten per cent. If enough phase space is available, 
the latter particles could then decay into the lighter chargino and neutralinos, Xii Xi an d 
X®, and neutral h,A,H or charged H ± bosons, with branching ratios again of the order 
of a few times ten per cent, making a total of a few percent branching ratio for the whole 
cascade. If the mass splitting between the lightest chargino or next-to-lightest neutralino 
and the LSP neutralino is large enough, another source of Higgs bosons will appear in a 
more direct way from the decays of xt ? X2 ~^ X?+ Higgs bosons, which has a rather large 
branching ratio. In some situations both possibilities are at work, leading to a substantial 
number of Higgs particles in the final state. [In the case of a relatively light charged 
Higgs boson, M#± < m t , there is an additional production mechanism, t — > bH + , with 
the t-quarks produced either directly or from the cascade decays of squarks and gluinos.] 
A key point in this analysis, is that the coupling of the Higgs bosons to chargino 
and neutralino states is maximal for higgsino-gaugino mixed states [16], while the gauge 
boson couplings to neutralinos are maximal for higgsino-like states. In the gaugino-like 
[i.e. when the higgsino mass parameter \n\ is much larger than the gaugino mass parameter 
M 2 ] or higgsino-like [i.e. in the opposite situation <C M 2 ] regions, this results into the 
dominance of the decays of the heavier charginos and neutralinos into the lighter ones 
and Higgs bosons, over the same decays with gauge boson final states in general. This 
is also the case of the little cascades in the gaugino region; the branching ratios for the 
decays of X2 an d Xi i nto the LSP and Higgs bosons, when kinematically accessible, are 
in general more important than for gauge boson final states. 

We will closely follow the analysis of Ref. [14], where the cross sections for squark 
and gluino production at the LHC have been discussed and where their decay branching 
ratios into the heavier charginos and neutralinos, which then decay into the lighter ones 
and charged Higgs bosons, have been analysed in detail. The analytical expressions of the 
various partial decay widths, including also the decays into neutral Higgs particles and the 
little cascades, have been given there. Here, we will simply show the total production cross 
sections times the decay branching ratios for the final states involving a single h, H, A and 
H ± particle in various scenarios. Throughout the analysis, we will assume the universality 
of the soft-SUSY-breaking gaugino mass parameters at the high-energy scale, leading to 
the relation m g ~ M 3 ~ 3M 2 ~ 6 Mi at low energies. 
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2.1 The various scenarios 

We will discuss the following four scenarios, which cover most of the possibilities for the 
cascade decays that we will discuss and span over the typical MSSM parameter space. 

— Scenario Scl: Here, the scalar partners of light quarks are heavier than gluinos 
[and top squarks]. For illustration, we chose rriq = 1.2m g = 720 GeV with the gluino 
mass taken to be m g = 600 GeV. The first- and second-generation squarks will then 
decay dominantly into quarks and gluinos, q — > gq [a smaller fraction will decay into 
quarks and chargino or neutralino final states]. In the case of lighter top squarks, since 
their mass is smaller than m g + m t , they will directly decay into chargino+bottom and 
neutralinos+top final states. The gluinos, which are directly produced or which come from 
the decays of squarks, will mainly undergo three-body decays into qqy® and qq'xt- This 
decay is dominantly mediated by the exchange of the lightest top squarks, which have a 
smaller virtuality [i.e. a mass smaller than the common scalar quark mass], thanks to the 
larger mixing between the stop eigenstates, which is proportional to m t . In this scenario, 
we will vary the higgsino mass parameter \x while the bino and wino mass parameters M\ 
and M 2 are fixed to M 2 = 200 GeV and M x = 100 GeV by our choice of m g , by virtue of 
the assumption of gaugino mass unification. 

— Scenario Sc2: This scenario is similar to the previous one, i.e. we assume that 
squarks are heavier than gluinos with the same ratio of trig = 1.2m g , but there is a major 
difference: here we take the gluino mass to be larger, m g = 900 GeV, which leads to 
M 2 = 2Mi = 300 GeV. For large enough \x values, the lighter charginos and neutralinos 
are gaugino-like with masses m x ± ~ m x o ~ 2m x o ~ M 2 . Thus for Ma ^ 130 GeV, the 
Higgs bosons are light enough, to render the mass splittings among the lighter gaugino 
states such that the decays xt ~^ X\H ± and X2 ~^ xV 1 ! H an d A can occur. Therefore, 
both the big cascade, which is the only one present in Scl, and the small cascade of eq. (2) 
are at work together in this scenario. 

— Scenario Sc3: Here gluinos are heavier than squarks, and we choose the common 
squark mass to be rriq = 800 GeV. Therefore, the decays g — > qq occur 100% of the time. 
Virtually, the electroweak cascades start with only the squark states, those coming from 
direct production and those from gluino decays. Here, we fix /i = 150 GeV and hence we 
are in the higgsino-like region, i.e. with a small value of the \i parameter with respect 
to M 2 [again universal gaugino masses are assumed at the high scale]. In this case, all 
squarks [in particular those of the first two generations whose couplings to the higgsino- 
like Xi)X\iX2 states are proportional to the small mass of their partner quarks and are 
therefore tiny] will mainly decay into the heavier chargino and neutralinos [which are 
gaugino-like with masses m x + ~ m x o ~ 2m x o ~ M 2 ]. For large enough M 2 values, there 
is sufficient phase space for the decay of the heavier gauginos into the lighter higgsino 
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states, with masses m x + ~ m x o ~ m x o ~ and Higgs particles with masses M$ < 200 
GeV to occur. 

- Scenario Sc4: This scenario is the same as the previous one as far as the squark 
and gluino sectors are concerned, but it is different for the electroweak gaugino sector. 
Indeed, here we choose a large value for the higgsino mass parameter, fx = 1 TeV, which 
now makes the lighter chargino and neutralinos gaugino-like. Squarks and gluinos will 
then decay into these states and, for large enough values of M 2 , the lightest chargino and 
the next-to-lightest neutralino can decay into the LSP and Higgs bosons [contrary to Sc3, 
where the little cascade was absent since xt,Xi2 where higgsino-like and degenerate in 
mass]. Of course, in this scenario, the big cascade is in principle kinematically possible 
but the branching ratios of first- and second-generation squark decays into higgsino states 
are strongly suppressed because of the small squark-quark-higgsino coupling. 

The Higgs sector will be treated using the program HDECAY, version 2.0 [17], which 
includes the important radiative corrections. The masses of the neutral h, H and charged 
H ± particles are obtained once the two input parameters at the tree-level, Ma and tan f3 
[as well as the scalar masses and other soft-SUSY-breaking parameters fi, Ah and A t , 
which enter the radiative corrections] are fixed. We will choose for illustration the input 
values: M A = 100,130 and 200 GeV and tan/3 = 10,30, which leads to the CP-even 
and charged Higgs boson masses given in Table 1. We will also discuss the scenario with 
Ma = 150 GeV and tan [3 — 5, where at the LHC only the h particle can be discovered in 
the direct production mechanisms, as discussed previously; the masses of the other Higgs 
bosons are also given in Table 1. 

To evade the LEP2 bounds on the light h boson mass 2 for small values of Ma and/or 
tan (3, we use a large value for the trilinear couplings A t > rriq, which increases the value 
of M h [the so-called typical-mixing scenario]. This leads to a large splitting in the stop 
sector: for rrig = 960 GeV and A t = 1.5 TeV for instance, one has ~ 840 GeV and 
mi 2 ~ 1.1 TeV [with a small variation with \i and tan/5], which also favours the cascade 
decays in Scl and Sc2, since the virtuality of the lighter t\ is smaller than that of the 
other squarks, enhancing the three-body decays of gluinos into higgsino final states. 

The produced neutral Higgs bosons will decay mainly into bb and t + t~ pairs, with 
branching ratios of respectively ~ 90% and ~ 10%, except for the CP-even h and H 
bosons when they are almost SM-like [for high values of tan f3 this occurs already when 
Ma ^ ( ^ )130 GeV for the h (H) states], where additional channels [such as WW^ 
decays] occur and will slightly suppress these rates. Depending on whether the tb threshold 

2 The absolute experimental bound from LEP2 searches of the h boson mass is Mh > 92 GeV, unless 
the Higgs boson is SM-like; in this case, the bound becomes Mh > 114 GeV [5]. This rules out small 
values of the parameter tan j3 in the absence of mixing in the stop sector. 
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M A [GeV] 


tan P 


M fe [GeV] 


M ff [GeV] 


M#± [GeV] 


100 


10 


97 


125 


127 




30 


99 


123 


126 


1 on 

130 


10 


116 


136 


151 






1 91 
1Z 1 


lOU 


101 


200 


10 


120 


201 


215 




30 


122 


199 


214 


150 


5 


111 


160 


169 



Table 1: Masses of the neutral and charged Higgs bosons for different inputs for Ma and 
tanp frrig = 960 GeV, A t = 1.5 TeV, /i = 100 GeV and M 2 = 270 GeV]. 

is reached or not, the H + bosons decay mainly into tb or t + u final states. 

The cross sections for squark and gluino production are evaluated as in Ref. [14], using 
the CTEQ3L [18] parametrization of the parton densities and with the scale defined as the 
average of the masses of the final sparticles. To be conservative, we will neither include 
the -fT-factors, which enhance the production cross sections [19], nor the possibly large 
QCD corrections to squark and gluino decays [20]. For the decay branching ratios, we use 
the average one defined in Ref. [14] for squark decays when they are heavier than gluinos 
[Scl and Sc2], i.e. that we sum over all possibilities for decaying left- and right-handed 
squarks as well as for up- and down-type squarks, keeping track of flavour and chirality 
to obtain a specific final state. For stops, we take into account the direct decays of the 
heavier ones into the lighter ones and Higgs or gauge bosons, but not the higher order 
decay modes [21]. For the three-body decays of gluinos, we take all possible channels 
into account and include masses for the third-generation fermions and sfermions and full 
sfermion mixing [22]. 

For charged Higgs bosons lighter than the top quark, we will also discuss their pro- 
duction from top decays. For moderate values of tan /3, 5 < tan j3 < 20, the branching 
fraction of the decay t — > bH + is rather tiny and the number of if ± bosons from the SUSY 
cascade decays is small. For larger tan (3 values, this branching fraction can be rather large 
[a few times ten per cent] and would lead to a substantial number of H^ 1 bosons from 
SUSY cascades. The top quarks are either produced directly, pp — > gg/qq — > tt or come 
from cascade decays of squarks and gluinos: two-body, g — > ttl 2 ,bbl <2 , or three-body, 
g — > ttXi,tbxti decays of gluinos and two-body decays of stops and sbottoms in top 
quarks and x states, U — > tx® and hi — > txj ■ 
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2.2 Cascades involving gauginos 

The variation of cross sections times branching ratios to obtain at least one neutral or 
charged MSSM Higgs boson in the final state from the big or little cascade, or from both, 
is shown as a function of fx in Figs. 1 and 2 for Scl and Sc2, respectively, and as a function 
of M 2 in Figs. 3 and 4 for Sc3 and Sc4. For each set, we take M A = 100, 130 and 200 
GeV for the top, middle and bottom rows respectively, and tan f3 — 10 and 30 for left and 
right panels, respectively. Fig. 5 shows the rates for a specific choice of Ma = 150 GeV 
and tan f3 = 5 in the above four scenarios. 

Let us first make a few general remarks on these figures. In all of them, we see that 
(jxBR for the four Higgs bosons can be rather large, exceeding the level of 0.1 pb in very 
large portions of the parameter space, and even reaching the level of ~ 10 pb in some 
cases. This means that more than 3 x 10 4 events, and up to 3 x 10 6 events, can be collected 
at the LHC for an integrated luminosity of / C ~ 300 fb -1 . In the cases of A, H and 
if ± bosons, the production rates are larger for smaller values of the masses [as expected 
from simple phase-space considerations], while they are stable for h boson production 
[since the variation of is mild for Ma =100-200 GeV]. The dependence of the rates 
on the value of tan j3 is not very pronounced in general, in contrast to the very strong 
tan /3-dependence for Higgs production in standard-like processes such as gluon-gluon 
fusion or associated production with bb pairs. 

In Scl with M 3 = 3M 2 = 6M 1 = 600 GeV, the mass differences between the lighter 
xt/X2 states and the Xi LSP are less than the minimum value of the lightest h boson mass 
[and thus also the A, H and if ± boson masses] considered, M h > 100 GeV, and therefore 
the little cascades, x° ~^ h, H, A + Xi and Xi ~^ + X\, are kinematically not allowed. 
The variation of <rxBR with fx is simply due to the variation of the branching ratios 
BR(<? -> xl,m, xf M') and BR(x3 i4 ,X2 = -»■ Xi,Xi,2 + Hi ggs) ; since the cross sections for 
squark and gluino production are constant, rrig and rriq being fixed. The production rates 
are larger for small (or moderate) \x values, \x < 400 GeV, where the gaugino (higgsino)- 
like states xf and Xs a have enough phase space to decay into the higgsino (gaugino)-like 
states xt and Xi,2- For larger values of fx, BR(g — > xt,^' xfqq') become smaller because 
of phase-space suppression, and crxBR drops out. 

For Ma ^ 150 GeV, the signal drops sharply for the CP-even and charged Higgs 
bosons for intermediate values of fx ~ M 2 ~ 200-250 GeV. In this case, the charginos and 
neutralinos are mixtures of gauginos and higgsinos and their mass differences are rather 
small, leading to phase-space-suppressed decays of the heavier states into the lighter ones 
and Higgs bosons. This is more pronounced for H^ 1 [and also H], which is the heaviest 
Higgs boson. However, the production rate for the CP-odd A boson is enhanced for small 
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M A values. This is due to a conjunction of several facts: the A boson is, together with h, 
the lightest Higgs particle and is therefore more favoured by phase space 3 ; its couplings 
to charginos and neutralinos are stronger than those of the h boson. Since the decays to 
Higgs particles share the branching ratio [together with the decays into gauge bosons], 
the suppression or the absence of decays into other heavier Higgs particles leads to an 
effective enhancement of the branching ratio for the decays into the pseudoscalar A boson. 

Note that for larger Higgs mass values, Ma ^ 150 GeV [say, Ma = 200 GeV], there 
is not enough phase space for decays into the heavier H, A, bosons in the small and 
intermediate /i range, fi < 250, GeV and only decays into the lighter h boson are allowed. 
For larger fi values, ax BR follows the same trend for H, A and if ± production, and even 
for h production for very large /z values [for H and A, this is expected since we are almost 
in the decoupling regime where they have the same masses and couplings]. The kinks 
followed by humps in the higher /i side are due to the opening of new decays channels, in 
particular channels involving the neutralino xh which is lighter than xf an d X% 

In Sc2, the cross sections times branching ratios for the four Higgs particles are smaller 
than in the previous case for low to medium values of the fi parameter, \i < 300 GeV, 
where only the big cascades are kinematically allowed. This is due to the larger gluino and 
squark masses considered, rriq = 1.2m g = 1080 GeV with m g = 900 GeV, which make the 
production cross sections of strongly interacting particles smaller. Nevertheless, ax BR 
is still rather large, being in general between 0.1 pb and a few picobarns for Ma ^150 
GeV [except in the case of H production, because of the smaller couplings to chargino 
and neutralino states] and the trend follows the one of the previous scenario. 

For larger values of the higgsino mass parameter, fi > 300 GeV, and since in this 
scenario M 2 = 2M\ = m g /3 = 300 GeV, we gradually enter the gaugino region where 
m x ± = m x o = 2m x o. Hence, there is now enough phase space for the lightest chargino 
and the next-to-lightest neutralino to decay into the LSP and a Higgs boson, at least 
for small Ma values. For Ma = 100 GeV, the little cascade occurs for all four Higgs 
bosons, but ax BR is much larger for h and H ± production [where it reaches the level 
of a few picobarns] than for the production of A and H, in particular for tan/3 ~ 10. 
This is mainly due to the strength of the chargino-neutralino couplings to Higgs bosons 
as discussed previously. For M A — 130 GeV, only the little cascades into neutral Higgs 
bosons are kinematically allowed, while for Ma ^ 150, only the little cascade into the h 
boson is kinematically allowed, the other Higgs particles being too heavy. 

3 In particular, the decay channel x° — f X?+Higgs can be kinematically open for the h and A bosons 
and closed for the H and if ± particles. Note also that the phase-space suppression is different for decays 
of charginos and neutralinos into CP-even and CP-odd Higgs bosons. 
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In Sc3, the cross sections times branching ratios for the production of the four Higgs 
particles from SUSY cascades are displayed as functions of M 2 , with the higgsino param- 
eter fixed to pi — 150 GeV. The common squark mass is chosen to be rriq = 800 GeV 
while the gluino mass varies with the gaugino mass parameter as rrig = 3M 2 . The values 
of (jxBR are rather large, being between 0.1 pb and a few picobarns, as soon as M 2 
exceeds the level of ~ 250 GeV. In this case the decay channels of the heavier chargino 
and neutralinos [which are gaugino-like with masses m x ± ~ m x o ~ M 2 ], into the lighter 
higgsino-like ones [with masses m x ± ~ m x\ 2 ~ ML are open only when |M 2 — \x\ > M$. 
This leads to a very large sample of signal events, ~ 10 5 with the expected luminosity 
of 300 fb _1 at the LHC. For larger M 2 values, the decay of the neutralino X3> with a 
mass m x » ~ Mi ~ \M 2 also opens up. For lower M 2 values, the mass splitting between 
the heavier and lighter chargino and neutralino states is too small, and the decays into 
Higgs bosons are not kinematically possible. Note that because the lighter chargino and 
neutralinos are almost degenerate in mass, the little cascades are absent. 

The production rates are larger for the h and H ± bosons: for h it is in general due 
to its lightness and consequently to having a more favourable phase-space, while for 
it is due to the fact that there are more possibilities for charged decays to occur 
[Xg 4 — > xfH T ,X2 ~^ XipH^ an d that charged current couplings are stronger than 
neutral current couplings. The production rates are relatively smaller for H and A boson 
production and we notice the equal production rate for large Ma values, when we are close 
to the decoupling limit. Note also that the production rates decrease with increasing M 2 
since, according to our assumption of universal gaugino masses at the high scale, trig 
increases and the cross section for gluino production drops out in the first place. 

In Sc4, the higgsino mass parameter is rather large: ii = 1 TeV. This leads to very 
heavy Xz 4 and xf [higgsino] states. The big cascades are disfavoured in this region for 
two reasons. First, squarks being relatively lighter [rriq = 800 GeV], they cannot decay 
into these heavy higgsinos. Secondly, three-body decays of heavy gluinos into the heavier 
higgsinos are suppressed by an extra factor of electroweak coupling-squared [compared to 
a two-body decay] and the gluinos mainly decay into squarks and quarks. Only the little 
cascades x 2 — ► h, H, A + x± and Xi H ± + x\ are possible if the Higgs bosons are light 
enough [the cascades are always possible for the production of the lighter h boson] . When 
these cascades occur, the cross sections times branching ratios are rather large, exceeding 
the level of 0.1 pb for not too large M 2 values and reaching the level of ~ 10 pb for light 
Higgs bosons and small M 2 values. For larger values of the gaugino mass parameter, 
cxBR drops following the decreasing phase space in the decay q — > x 2 <?, Xi <?'> while a 
further dampening is caused by the lower gluino production rate due to a simultaneous 
increase in gluino mass. 
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Figure 1: Cross sections times branching ratios for Higgs production in Scl. 
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Figure 2: Cross sections times branching ratios for Higgs production in Sc2. 
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Finally, in Fig. 5, we show the cross sections times branching ratios in the four selected 
scenarios in the case Ma = 150 GeV and tan/3 = 5. Since, as we discussed earlier, the 
dependence on the value of tan (3 is not very pronounced, the trend is half-way between 
what occurs for Ma = 130 GeV and for Ma = 200 GeV. Except for the regions of 
intermediate /i value in Scl and Sc2 and small M2 values in Sc4, ax BR for the production 
of the heavier Higgs bosons H, A and H ± via the cascades are substantial, exceeding the 
level of 0.1 fb in large regions and even reaching the level of a few picobarns in some cases. 




300 350 400 450 500 300 350 400 450 500 



Figure 5: Cross sections times branching ratios for MSSM Higgs boson production for 
Ma = 150 GeV and tan(3 = 5 in the four scenarios Scl to Sc4- 
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2.3 Higgs bosons from direct decays of stops and sbottoms 

If the mass splitting between two squarks of the same generation is large enough, as is 
generally the case of the (t,b) iso-doublet, the heavier squark can decay into a lighter 
one plus a Higgs boson $ = h, H, A [or a Z boson]. If, in addition, there is enough mass 
splitting between the stops and sbottoms, it is possible for the heavier one to decay into 
the lighter one and H ± [or W^] states. The case of charged Higgs bosons has been studied 
in detail in Ref. [14] and we concentrate here on the stop and sbottom decays into the 
neutral Higgs particles, t 2 — > fi$ and b 2 — > &i<E>. 

While the production mechanisms for stops and sbottoms are mainly governed by QCD 
and the cross sections depend only on the mass of these particles [as well as on the gluino 
mass when the decays g — > tt, bb are kinematically possible], their subsequent decay rates 
into Higgs particles depend on many parameters. First, since these squarks are heavy, 
they have many possible decay modes besides the ones into Higgs bosons discussed here: 
decays into charginos/neutralinos and quarks, which in general have rather sizable decay 
widths, as well as decays into gauge bosons. In addition, these decays depend strongly 
on the Higgs-squark couplings, which in the case of the neutral Higgs particles are given 
by [see Ref. [14] for a discussion of the couplings and decay widths]: 

g qiq2 h oc M|sin2^(2/ 3 -4e g s^)sin(a + /3) + 2m g cos2^(A (? ^ + 2/ 3 /i^) 



9q\qiA — 9q-2,q\A ^ m q 



9 qi q 2 H oc -M\ sin 20,(2/J - 4e,s^) cos(a + P) + 2m q cos 26 q {A q r\ - 2/J /1 r q 2 ) 



pi + A q (tan f3)~ 



(5) 



where a is the mixing angle in the CP-even Higgs sector of the MSSM, 9 q is the squark 
mixing angle [for 9 q = 0, qi = q~R in general], and the coefficients r\ 2 are the normalization 
factors of the h, H boson couplings to fermions normalised to the SM Higgs boson coupling. 
The largest component of the CP-even Higgs boson couplings is the one proportional to 
the quark mass [in the case of sbottom couplings, a tan (3 enhancement is present for the 
coupling of the H state]. However, if the mass splitting between the two squarks is large, 
the squark mixing angle approaches the value \6 q \ ~ 7r/4 and the term cos20 g approaches 
zero [the other component of the coupling, the one proportional to M z , although maximal 
since sin20q ~ 1, is not large enough to make the Higgs-squark coupling strong]. In the 
case of no mixing, where the coupling becomes strong, there is in general not large enough 
mass splitting between the two squarks [at least if one assumes common scalar mass terms] 
to make the decays into Higgs bosons possible. Thus, the only decays that could have 
sizeable partial widths would be the channels t 2 (b 2 ) — > ti{pi)A whereas, in the case of 
sbottoms, a large value of tan/5 is needed to make the sbottom mass splitting large for 
the decay to occur and to compensate for the smallness of m h in the coupling to have a 
sizeable rate. 
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Figure 6: Cross sections times branching ratios for the production of the pseudoscalar 
Higgs boson from direct decays of heavier stop and sbottom squarks, £2(^2) — ► A + tiipi), 
for the values Ma = 130 GeV and tanfi = 10 (left panel) and 30 (right panel). The full 
lines correspond to fa, 62 squarks originating from decays of gluinos produced in pp — > qg 
and gg processes, while the dashed lines correspond to top and bottom squarks produced 
directly in the processes pp — > ^2^2- 
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In Fig. 6, we present the cross sections times branching ratios for the production of 
a single A boson in these cascade decays 4 . We have taken into account all production 
channels for gluinos and squarks, including the third generation and all possible decays of 
squarks, i.e. decays into \ states and quarks, decays into gauge bosons and the possible 
decays into H ± bosons. Representative values of the parameters are: rriq = 600 GeV 
assuming all soft squark masses are degenerate, m g = 3M 2 , A q = 1.5 TeV and tan/3 = 
10, 30; we fixed the pseudoscalar Higgs mass to Ma = 130 GeV and varied M 2 or /x. 

The upper figures show erxBR as a function of /x for the fixed value M 2 = 400 GeV. 
The solid lines represent the contribution from g cascade decays with the production 
processes gg and qg, whereas the dashed lines are for direct production of q 2 q 2 pairs [note 
that for tan f3 = 10, the sbottom contribution is absent as the mass splitting between b 2 
and bi is not sufficient to access the decay b 2 - ► b\ +A}. The rates increase with increasing 
/x values for two reasons: first, the decay rates t 2 — > t\ + A increase with /x, as can be seen 
from the Atit 2 coupling discussed above; secondly, a smaller phase space is accessible to 
stop decays into the heavier charginos and neutralinos, at higher values of /x. [Note that 
the production cross sections are not sensitive to /x, except for a mild dependence of the 
t 2 mass on this parameter.] For the higher tan (3 case, the situation is similar, except that 
at a certain value of /x (~ 700 GeV), the splitting between m~ b2 and rag becomes large 
enough to permit the opening of the channel b 2 — > b\ + A, which eventually contributes 
in both the g cascade decay and in the direct production of b 2 . 

The two lower figures show crxBR as functions of M 2 with a fixed value of /x = 
1000 GeV. In the low tan (3 case, the decay of b 2 — > b±A is not accessible [as we have seen 
earlier], whereas for the high tan (3 case it contributes at the same level as t 2 . Here, the 
rate out of the g cascade decay falls very sharply as the gluino production cross sections 
decrease with increasing M 2 = \rrig. On the other hand, the stop and sbottom production 
cross sections are almost independent of M 2 ; their slightly rising behaviour is due to the 
increased decay branching ratios in the A boson because of the suppressed rates in the 
gaugino decay channels, which are more phase-space-disfavoured for larger M 2 . 

Thus, the production of A bosons from direct decays of t, b squarks [as well as the 
production of bosons from the decays b — > t\H~, which were discussed in Ref. [14] and 
found to lead to similar results] , although in general with much smaller cross sections times 
branching ratios than Higgs production from cascades involving charginos and neutralinos, 
might lead to the detection of these particles in favourable situations. An interesting 
feature is that these processes could allow, in principle, to disentangle the H and A bosons 
in the decoupling limit, which is otherwise very difficult to achieve in other processes. 

4 We explicitly verified that, for the CP-even neutral Higgs boson channels, the rates are suppressed 
by few orders of magnitude and are too small to be detected at the current luminosity options at LHC. 
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2.4 Charged Higgs bosons from top decays under SUSY cascade 

The decays of the top quarks that are pair-produced at the LHC are known to be a 
very important source of charged Higgs bosons for m t > m H ± . Situations permitting, as 
we will see in this section, cascades of SUSY particles leading to top quarks followed by 
t — > bH^ may also become a significant source of iJ ± bosons. 

Charged Higgs bosons can be produced in SUSY cascade decays via the pair production 
of gluinos, top and bottom squarks, as in eq. (4), followed by their cascades to top quarks, 
which subsequently decay to charged Higgs bosons. To have a broad estimate of how 
significant this new source of charged Higgs particles can be, we discuss cases with the 
four reference scenarios Scl-Sc4 presented previously, and calculate the rates for the 
production of a single charged Higgs boson final state from these cascades. To simplify 
our discussion, we will not consider the possibility of production via the bg production 
mechanism [which has a very low cross section] as well as from decays of heavier chargino 
and neutralino states, x ~ * t + X [which would simply increase the sample]. 

In Fig. 7, we present the results of our study for the four scenarios with the choice 
tan/5 = 30 and Mu± = 130 GeV, in terms of variations of cross sections for different 
production processes, times the branching ratios for the decays g,q — > t+X and t — > bH^. 
The variations are presented with respect to /i for Scl and Sc2 and with respect to M2 for 
Sc3 and Sc4. We do not present the results for a lower tan (3 value since it can be estimated 
in a straightforward manner from the one we are analysing: while the production cross 
sections and decay branching ratios of SUSY particles remain almost constant, the final 
ax BR drops by a factor ~ 1/ tan 2 /? [i.e. by a factor ~ 9 when going from tan/5 = 30 
to 10] because of the suppressed Yukawa coupling, gn±tb ^ ^btan/5 [the component of 
this coupling that is proportional to m t is suppressed by a factor 1/ tan/?]. Therefore, the 
total ax BR in the tan/5 = 10 case is one order of magnitude smaller than in Fig. 7, while 
the relative magnitude of the various contributions is approximately the same. Note that, 
in each figure, the long arrow along the y axis represents the H ± production cross section 
via direct tt production. 

The axBRs are suppressed by roughly an order of magnitude when going from Scl 
to Sc2 owing to an increased gluino mass and consequently larger squark masses. The 
variations with respect to /i present in these scenarios are mainly determined by the 
resulting variation of the tt\\ and btxf couplings. For small /x values [«C M 2 ] the lighter 
neutralinos and chargino have masses m x o^m x ± ~ /1 while they are higgsino-like, yielding 
larger branching ratios in those g decay channels and also for t\ — > tx°j{j =1-4). However, 
at a certain value of /i, xl become so heavy that they are not accessible in stop decay; 
this results in a smaller or vanishing BR(ti — > ^34) value. A further increment of crxBR 
beyond this /1 value in Scl is due to a gradual increment of the gaugino content of Xi,2- 
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The profile of the g contribution in Scl roughly follows that of tit *, because the same 
decay mode ti — > tx± is involved there, to end up with top quarks. The flatness of ax BR 
for b\b \ production then indicates the low sensitivity of the bx ± t coupling to //, which is 
due to the small value. For Sc2 the trend is similar, except for the fact that now, 
the kinks do not appear in the frame. This is due to an increased value of m g [and 
hence, of rriq = 1.2rrig] which effectively pushes to higher values, enabling access to 
all neutralino states in its decay t\ — > tx°j for the entire range of /x. 

In Sc3 and Sc4 the contributions from gg and gq production drop naturally with 
increasing M 2 = \M gi as a result of the smaller production cross section for heavier 
gluinos. The peaks at about M 2 ~ 270 GeV indicate the onset of the two-body decay of 
a gluino, g — > tt\. While going from Sc3 to Sc4, we see that the situation does not change 
much, except for a slight increase of ax BR in the t\ pair-production case and a drop in 
the b\ pair-production case beyond 500 GeV in Sc4. The flatness of the t\t\ contributions 
in Sc3 indicates that the main contribution is coming from t\ — > tx® where % —1 and 2, 
i.e. from the lighter neutralinos with masses close to /i. In the case of sbottoms, the main 
decay mode seems also to be bi — > txi , this chargino being a higgsino. 

With /i — 1 TeV for Sc4, the decay bi — > txi gradually drops with increasing M 2 , 
leading to a decrease of available phase-space [here, the lighter chargino is gaugino-like] 
and gets closed completely for a high M 2 value. However, in the stop case there is an 
increase in ax BR at about 500 GeV for these contributions, which indicates the onset 
of a gradual increment of the branching ratio t\ — > txi, a ^ the expense of other channels 
including decays into higher neutralinos and charginos as they are getting closer to phase- 
space suppression with the increase of M 2 for a value of fi as high as 1 TeV. 

Thus, for the large value of tan (3 = 30 considered and for Scl and Sc2, this additional 
contribution to single charged Higgs production may even become comparable to the one 
from the cascades involving charginos and neutralinos, in a rather conservative estimate 
[the cross sections times branching ratios being at the level of 0.1 pb]. For Sc3 and 
Sc4, these new contributions can only be significant for low gluino masses and naturally 
drop sharply with the increase in the latter. As mentioned already, and contrary to the 
processes discussed in section 2.2, the event rates decrease quadratically with tan/3 and, 
for instance, are an order of magnitude smaller for tan f3 = 10. 

In short, for not too heavy a gluino, m g < 1 TeV, and large values of tan tan (3 > 20, 
top quarks coming from SUSY particle cascade decays may contribute heavily to charged 
Higgs boson final states in some regions of the MSSM parameter space. This would always 
strengthen the signal, although one has to be rather meticulous about the backgrounds 
[in particular, it would not be an easy task to separate these events from the ones where 
the H ± bosons originate from the top quarks, which are directly produced]. 
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Figure 7: Cross sections times branching ratios for charged Higgs boson production from 
decays of top quarks originating from cascade decays of strongly interacting SUSY particles 
in the scenarios Scl to Sc4 and for the values Ma = 130 GeV and tan/3 = 30. For 
tanj3 = 10, the trend is similar except that the rates are one order of magnitude smaller. 
The arrows on top of the figures correspond to the cross section for gg/qq —* tt times the 
branching ratio for the decays t — > H + b and i H~b. 
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3. Detector simulation analysis 

In the previous section, the production mechanism of MSSM Higgs bosons through 
strongly interacting SUSY particle cascade decays has been extensively discussed for sev- 
eral SUSY scenarios. The important feature of this mechanism is the quasi-independence 
of the production cross sections times branching ratios on the value of tan/5. This is in 
contrast to the strong tan j3 dependence of the direct production mechanisms for the heav- 
ier charged and neutral MSSM Higgs bosons, which in the standard processes [gg — > H/A 
or gg/qq — > bbH/A for the neutral Higgs bosons and gg — > t&iJ or gb — > H~t for the 
charged Higgs particles] lead to a significant number of events only for larger value of this 
parameter, for which the Hbb, Abb and H^tb couplings are sufficiently enhanced. On the 
other hand, the masses of the MSSM Higgs bosons accessible in the SUSY cascade de- 
cays will probably be limited to the range Ma ^ 250 GeV by phase-space considerations. 
Therefore, one might hope to cover the lacuna in the discovery reach for the heavy MSSM 
Higgs bosons at the LHC, which is left after including the SM [and possibly sparticle] 
decay modes of these particles, i.e. for low pseudoscalar Higgs boson masses, 120 GeV 
< Ma ^ 200 GeV and moderately small tan/3 values, 3 < tan/3 < 10. 

In this section, we will present a first analysis that will demonstrate the detection 
potential of the MSSM Higgs bosons in SUSY particle cascade decays. As already men- 
tioned, the subject is rather involved and in the present paper, we will simply perform 
a fast Monte-Carlo simulation of the production signal and the various SM and SUSY 
backgrounds for H, A and H ± boson production via squark and gluino decays through 
the big and small cascades. [We will not discuss direct decays of top and bottom squarks 
into Higgs particles or the production of H ± bosons through top quark decays]. For il- 
lustration, we will focus on the case of tan/3 = 5 and Ma = 150 GeV, which leads to 
a light Higgs boson with a mass of M h ~ 110 GeV [since h is not yet SM-like for this 
set of parameters, this mass value is not yet excluded by LEP2 searches] and a charged 
Higgs boson with a mass M H ± ~ 170 GeV [a value for which the decays t — > H + b would 
be strongly suppressed], in the scenarios discussed in section 2. The cross sections times 
branching ratios for these scenarios have been displayed in Fig. 5. 

The signal and the background events are generated with HERWIG 6.4 [26]. The SUSY 
spectrum is calculated 5 with ISASUSY 7.58 [27] and then interfaced to HERWIG using the 
ISAWIG package [28]. To be as realistic as possible, we also simulate some aspects of the 
CMS detector using CMS JET 4.801 [23], which contains a parametrized description of the 

5 Some of the sparticle and Higgs boson masses as well as some of the branching fractions for the 
cascade decays obtained with this program are slightly different from those obtained in our numerical 
discussion of the previous section. However, this discrepancy does not alter the main features of our 
analysis and can be resolved [i.e. one can obtain a similar spectrum] by slightly changing the inputs. 
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CMS detector response. The effects of event pile-up at the LHC have not been included, 
but are expected to be small for our analysis. The features that will allow the distinction 
between the Higgs boson signals from the SM and SUSY backgrounds will be discussed 
in the next two subsections, and we will argue that the SM backgrounds can be efficiently 
suppressed. However, a more difficult task will be to distinguish the signals from other 
SUSY cascade decay processes since all squark and gluino production at the LHC will 
end up in lightest neutralinos and fermions through gaugino and slepton decays. Hence, 
a good understanding of the nature of these SUSY backgrounds will be needed prior to 
any search for MSSM Higgs bosons in cascade decays. 

3.1 Neutral Higgs bosons 

We start by describing a first-analysis strategy that will allow the observation of the 
neutral h, H and A bosons produced in the squark/gluino cascades. Two decay modes of 
these neutral Higgs particles are promising, namely A,H^ bb, with a branching ratio 
around 90%, and A,H — > rr, with a branching ratio around 10% [the branching ratios are 
slightly smaller in the case of the lighter h boson]. Because of the dominant branching rate, 
we will now focus on the bb decays, although the r + r~ decay may have some advantages 
due to the lower multiplicity of r-jets with respect to 6-jets in the cascade background. 
Eventually both modes should be studied since the ratio of their signals could be used to 
support the evidence for a Higgs signal in the SUSY cascade decays. 

In all scenarios, the h, H and A bosons originate from the gg, qq, qq* or qg production 
processes through the chain illustrated in eq. (1) or (2). In Scl and Sc3, only the first 
chain is allowed, i.e. only the "big cascade" mechanisms are at work. In Sc2, the light h 
boson can be produced through the "little cascade" mechanism of eq. (2). In Sc4, also the 
H, A bosons can be produced in the little cascades. In all cases the background consists 
of the same gg, qq, qq*, qg processes, where the squarks and gluinos decay into neutralinos 
and charginos that do not lead to Higgs bosons in their subsequent decay chains. 

We will study the features of the signal and the backgrounds in order to determine a 
set of selection criteria that should serve the following generic purposes: 

• They should make sure that the events can pass the CMS on-line trigger levels with 
a high efficiency. 

• They should reduce the SM backgrounds to the level where they are negligible with 
respect to the SUSY backgrounds. 

These two purposes are closely related, since the triggers are designed to reject SM 
jet backgrounds. In contrast to SM processes, cascade decay Higgs events contain a large 
number of very energetic jets and a large transverse missing energy. By selecting only 
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events with these features using the criteria that we will propose, trigger efficiencies of 
more than 90% can be reached [24]. Additional cuts intending to suppress the SUSY 
background i.e. the SUSY events not containing any Higgs bosons, will depend on the 
details of the kinematics of this background. In the present paper, to be as unbiased 
as possible, we will not adapt any special cuts to reject this SUSY background. After 
suppression of the SM background, we will merely look for resonances in the bb invariant 
mass spectrum for the four scenarios, comparing this spectrum to its equivalent in the 
case where Higgs boson production in the SUSY cascades is not allowed. 

We will now discuss the kinematical distributions of the SUSY cascades and the SM 
background in order to obtain a recipe that allows us to suppress the latter while preserv- 
ing the former. With 'SUSY cascades' we mean the production of all sparticles according 
to cross sections obtained from the four scenarios as described in the previous section. 

The concrete parameter choices for these scenarios are [the gluino mass m g has been 
slightly changed in Sc3 and Sc4 to obtain a spectrum that is closer to the one discussed 
in section 2, while the common slepton mass has been fixed in all four scenarios to the 
value m z ~ = 500 GeV]: 

• Scl: M 2 = 2Mi = 200 GeV, fi = 300 GeV, m- g = 600 GeV and m~ q = 720 GeV. 

• Sc2: M 2 = 2Mi = 300 GeV, fj, = 450 GeV, m g = 900 GeV and m q = 1080 GeV. 

• Sc3: M 2 = 2Mi = 350 GeV, fx = 150 GeV, m g = 1200 GeV and m q = 800 GeV. 

• Sc4: M 2 = 2Mi = 350 GeV, fi = 1000 GeV, m g = 1200 GeV and m q = 800 GeV. 

As SM backgrounds, we have only simulated the main tt process. Also other QCD 
processes [bb, cc and light quark production] will contribute to the SM background; how- 
ever, it is very difficult to produce a reliable QCD sample in this case, since extreme 
kinematical fluctuations are necessary for this type of background to be within the signal 
cuts. Therefore we have not included the QCD background, being confident that any set 
of criteria that is able to suppress the tt background will also be effective against the other 
QCD background. We recall that ISASUSY 7.58 and HERWIG 6.4 are used for the event 
generation and that the detector aspects were simulated using CMS JET 4.801 [23], which 
contains fast parametrizations of the CMS detector response and a parametrized track 
reconstruction performance based on GEANT simulations for 5-quark tagging [25] . Results 
for signal and background are summarized in Figs. 8-11. 

In Fig. 8 the jet multiplicity is shown, for the SUSY cascades and the SM tt back- 
ground, for each of the four scenarios. The first two scenarios have a jet multiplicity 
distribution that peaks around 8 jets, while the two last scenarios show somewhat lower 
jet multiplicities. This is linked to the hierarchy of squark, gluino and gaugino masses as 
described in section 2. For the SM tt background, the distribution peaks around 5 jets. 
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Figure 8: Normalized distribution of the jet multiplicity in the SUSY cascades (full line) 
and the SM tt background (dashed line) events, for the four scenarios. Typically, the 
SUSY cascades events contain more jets than the SM background events. 
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Figure 9: Normalized distribution of the of the hardest jet in the SUSY cascades 
(full line) and the SM ti background (dashed line) events, for the four scenarios. The 
distributions peak around 100 GeV for the SM background and above 300 GeV for the 
SUSY cascades. 
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Figure 10: Normalized distribution of the transverse missing energy in the SUSY cascades 
(full line) and the SM ti background (dashed line) events, for the four scenarios. The 
distribution of the SM background peaks at low values and drops steeply, while the SUSY 
cascades has a very broad distribution reaching very high values. 
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Figure 11: Normalized distribution of the effective mass consisting of the sum of the of 
the jets and the E™ 1SS in the SUSY cascades (full line) and the SM tt background (dashed 
line) events, for the four scenarios. The distributions of the signal and the SM background 
are clearly separated in all scenarios. 
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Requiring a large number of jets in the event will enhance the sparticle production versus 
the SM background in most scenarios. 

In Fig. 9, we can see the transverse energy of the hardest jet in the event, for the 
SUSY cascades and the SM tt background events, in each of the four scenarios. The 
distribution peaks around 100 GeV for the SM background and above 300 GeV for the 
SUSY cascades. Demanding the E T of the hardest jet to be above ~ 300 GeV will 
therefore strongly suppress the SM background. 

Fig. 10 shows the transverse missing energy in the events, for the SUSY cascades 
and the SM tt background, in each of the four scenarios. The distribution of the SM 
background peaks at low values and drops steeply, while the SUSY cascades have a very 
broad distribution reaching very high values. Again, requiring a large E^ 1SS in the events 
will strongly favor the sparticle production with respect to the SM background. 

In Fig. 11, the variable E^ tal is shown, consisting of the sum of the Ex of all the jets 
and the E™ 1SS in the event. The distributions of the SM ti background and the SUSY 
cascades are clearly separated in each of the four scenarios. Selecting only events in which 
Efi tal is above ~ 1200 GeV will help suppress the SM background. This selection cut in 
itself is very effective, however it becomes somewhat redundant if applied after the two 
previous cuts on E^ and E™ 1SS . 

Combining the above observations allow us to strongly suppress the Standard Model 
background while preserving most of the SUSY cascade events. As a next step, we will 
study the bb invariant mass spectrum and look for resonances of the h (110 GeV), A (150 
GeV) and H (160 GeV) Higgs bosons. In order to perform this, we will demand that 
there be at least two jets in the event carrying a significant 5-tag and having E T values 
compatible with the ones originating from a Higgs boson. As mentioned before, no other 
cut will be made to suppress the SUSY background in view of the uncertain nature of 
this background. We can thus summarize the selection strategy as follows: 
i) We require the event to contain at least 5 jets. 

ii) The hardest jet in the event should have an energy E T larger than 300 GeV. 
Hi) The transverse missing energy E™ 1SS should be larger than 150 GeV. 
iv) The variable E$ tal = ££^ cts + E^ s should have a value larger than 1200 GeV. 
v) The event should contain at least two 6-jets with 45 < E T < 120 GeV. In order for 
the jets to be 6-tagged, we demand that both jets contain at least two tracks with 
a significance of the transverse impact parameter cr{ip) = ip xy /Aip xy > 3. 
If more than two jets in the event fulfill the last condition, the bb invariant mass is 
calculated using the two 6-jets that are closest to each other in r/-(f) space. Applying the 
above condition a{ip) > 3 leads to a b-tagging efficiency of about 40%, with a mistagging 
probability of less than 0.5%. 
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For each of the four scenarios, we will apply this selection strategy to two different 
squark/gluino samples: one where all decays are allowed (SUSY signal) and another one 
where the chargino/neutralino decays into Higgs bosons are switched off (SUSY back- 
ground). By overlapping the bb invariant mass distribution obtained with both samples, 
we can study the effects of Higgs bosons in the SUSY cascades in an unbiased way. The 
main results are summarized in Figs. 12-15, assuming a luminosity of / C = 30 fb _1 . 

Fig. 12 shows the bb invariant mass spectrum for the SUSY signal overlapped with 
the SUSY cascade background for Scl. Both signal and background feature a Z boson 
peak. In our fast simulation code, the "measured" M b i is typically underestimated by 
about 10% compared to the "true" M b g. However, the Z peak seems to be shifted to 
lower values by more that. This is due the large presence of x° ~~ > ^Xi decays in the 
cascades, the bb invariant mass distribution of which features a kinematical edge around 
80 GeV. This edge interferes with the Z peak and leads to the apparent peak around 
70 GeV. There are two more, relatively small, peaks visible in the signal M hb spectrum: 
one around 100 GeV and one around 140 GeV. The first one corresponds to the presence 
of the h boson while the second one corresponds to the A, H bosons. In this scenario, 
all Higgs bosons are being produced in the big cascades, i.e. they originate from heavy 
chargino/neutralino decays. In order to obtain a more convincing S/B ratio, this scenario 
requires additional selection cuts to suppress the SUSY background. In this figure, also 
the Standard Model tt background is shown. Clearly it is small compared to the SUSY 
background, justifying our motivation to neglect the QCD background. 

Fig. 13 shows the M b i distributions for Sc2 and one can see that a large signal peak is 
coming from h — > bb decays. This is due to the fact that the h boson is being produced 
in the little cascades [x° — ► hx° has a branching fraction of 96%], while the H and A are 
only produced in the big cascades. Therefore the H, A peak is less visible. Moreover, 
there is a broad smearing due to the large combinatorial background coming from the 
wrong pairing of 6-jets since the H or A state is often accompanied by a X21 which decays 
X2 —> hxi, leading to four b quarks in the event. In the background, only the Z boson 
peak can be seen. 

In Fig. 14, the result is shown for the third scenario. There, as can be seen in Fig. 3, 
the production rates for all Higgs bosons are rather high. There are no little cascades, and 
the rates for h and H, A production are therefore not very different. A double Gaussian 
fit clearly shows the corresponding peaks, while in the background only the Z boson peak 
is visible. 

In Fig. 15, the bb invariant mass spectrum is shown for Sc4. In this case, not only the 
lighter h boson but also the heavier H and A particles are produced in the little cascades, 
leading to very clear peaks for the light and the heavy neutral Higgs bosons. Again, in 
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the background spectrum, only the Z boson peak can be observed. Both this scenario 
and the previous one feature squarks lighter than the gluino [which has a mass around 
1 TeV] ; from the figures, it is clear that this condition strongly enhances the visibility of 
the neutral Higgs boson signals. 

Conservatively estimating the significance of the peaks as S = S/ VB, where S is the 
number of events in the peak and B the number of events "below" the peak [i.e. below a 
hand-drawn line between the two shoulders of the peak], we obtain a 5<r significance for 
the A, H peak in Sc4. With 30 fb _1 of integrated luminosity, enough data are collected to 
observe a 5<r excess also for Sc3. Scl and Sc2 will require a more detailed study and/or 
a larger integrated luminosity. 

This simple analysis of the neutral Higgs bosons produced in the SUSY cascades 
shows that this alternative production mechanism looks very promising. The Standard 
Model backgrounds can be efficiently suppressed by the selection criteria outlined above. 
By studying the bb invariant mass spectrum, we have shown that for all the considered 
scenarios the neutral Higgs bosons are visible after a few years of low luminosity running 
of the LHC collider. 

This has been exemplified in the case where Ma = 150 GeV and tan/5 = 5. However, 
as discussed in section 2, the dependence of the signal cross sections times branching 
ratios on the parameter tan (3 is rather weak in the cascade decays and we can therefore 
extrapolate this conclusion to all values of this parameter. For the Higgs mass range that 
can be probed in these processes, we are limited by two factors: 

(i) In the low mass range, M A > 90 GeV, we get rather close to the Z boson peak 
and the signal and background overlap; however, the cross section times branching ratios 
for the signal can be rather large in this phase-space-favoured case and the signal peak 
can be much larger than the background peak. Moreover, by normalizing the Z — > bb 
distribution to the Z — > ee, distribution, any excess can be clearly established. 

{%%) For larger Ma values, the A, H bosons become too heavy and phase space becomes 
penalizing [the branching ratios for the decays of heavier chargino/neutralino states into 
the lighter ones and the H, A bosons become too small]. However, in many scenarios one 
can still have reasonable Higgs production rates up to pseudoscalar masses of the order 
of M A = 250 GeV. Certainly, in Scl-Sc4 discussed here, one can probe masses up to 
Ma = 200 GeV, since there still are reasonable event rates, as shown in Figs. 1-4. 

Therefore, it appears that the heavier MSSM Higgs bosons originating from cascade 
decays of squarks and gluinos can be detected at the LHC for any value of tan f3 and for 
pseudoscalar Higgs boson masses approximately in the mass range between 100 and 200 
GeV [the lightest h boson can be detected in its entire mass range] in a few years of LHC 
running with a moderate luminosity. 
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Figure 12: Distribution of the bb invariant mass for the SUSY signal events on top of 
the SUSY cascade background assuming scenario Scl (with 30 fb^ 1 luminosity). At the 
bottom, the distribution for the Standard Model ti background is plotted. Two rather 
small peaks can be seen in the signal, corresponding to the presence of h and A, H bosons 
originating from the big SUSY cascades. 
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Figure 13: Distribution of the bb invariant mass for the SUSY signal events on top of 
the SUSY cascade background for scenario Sc2 (30 fb^ 1 ). The SM ti background is also 
shown. A large signal peak is visible, corresponding to the h bosons that are abundantly 
produced in the little cascades; also a much smaller and broader peak can be seen, signaling 
the presence of A and H bosons coming from the big cascades. 
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Figure 14: Distribution of the bb invariant mass for the SUSY signal events on top of 
the SUSY cascade background for scenario Sc3 (30 fb^ 1 ). The SM ti background is also 
shown. Two peaks can be distinguished, a large one corresponding to the h boson and a 
smaller one corresponding to the A, H bosons, all originating from the big cascades. 
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Figure 15: Distribution of the bb invariant mass for the SUSY signal events on top of 
the SUSY cascade background for scenario Sc4 (30 fb^ 1 ). The SM ti background is also 
shown. Two peaks are clearly visible, corresponding to the h and A, H Higgs bosons, all 
being produced in the little cascades. 
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3.2 Charged Higgs bosons 

After having demonstrated the potential for observing the neutral Higgs bosons in the 
SUSY cascade decays, we will now study the charged Higgs production through this 
mechanism. Since in all considered scenarios, Ma — 150 GeV, the charged Higgs boson 
mass is close to the top quark mass, M H ± ~ 170 GeV, and top quark decays into bottom 
quarks and H^ 1 bosons will be strongly suppressed by phase space. For the same reason, 
the charged Higgs decay mode to tb will be suppressed and the if ± particles will decay 
into tv final states with a branching ratio of 95% [some small fraction would also decay 
into cs final states for low values of tan (3} . This enforces our choice of the decay mode of 
investigation to be H ± — > 

We will select events containing a hard r-jet plus additional hard jets (often 6-jets) 
accompanied by a large amount of missing energy due to the presence of the lightest 
neutralinos and the neutrinos. Since there are numerous sources of missing energy in 
the cascade decays, an invariant mass reconstruction of the charged Higgs using the rv T 
mode is not possible. Also, in the tb decay mode, if the branching ratio had been sizeable, 
a reconstruction of the 4-jet system [2 6-jets + 2 jets from the W boson] would have 
suffered from a huge combinatorial background in the environment of the cascade decays. 

Though deprived from the powerful tool of invariant mass reconstruction, the hadronic 
H ± — > t ± v t mode has a special feature that can be exploited in distinguishing it from 
the main background decay W ± — > r ± z/ r : the tau polarization effects for both channels 
are quite different [30]. 

Indeed, due to the scalar nature of the decaying H + particle, the r + lepton from 
its decay is produced in a left-handed polarization state. In the simplest scenario of 
t + — > ti + u decay, the right-handed v is preferentially emitted in the direction opposite 
to the r + in the tau rest frame, so as to preserve the polarization. On the contrary, in 
the background, the r + from the W + — > t + v decay is produced in a right-handed state 
due to the vector nature of the W + boson, forcing the v from r+ — > tt + u to be emitted in 
the same direction as the r + in the tau rest frame. Therefore, harder pions are expected 
from the signal than from the W bosons in the background. By selecting events in which 
the fraction of the r-jet transverse momentum carried by the charged pion is large, the 
SM and SUSY background involving W bosons can be suppressed relative to the signal. 
The tau leptons in the background coming from either f and Z decays or the decays of 
the MSSM neutral Higgs bosons [in particular H and A] can however not be suppressed 
this way. 

Thus, as a signature for the charged Higgs boson in the SUSY cascade decays, we will 
look for highly energetic r-jets in which the leading charged particle carries most of the 
momentum. 
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For the simulations, we have again used ISASUSY 7 . 58 [27] for the calculation of the 
SUSY particle spectrum and the HERWIG 6.4 [26] Monte-Carlo event generator for the 
signal and background generation. This generator takes into account the spin correlations 
between the production and decay of all heavy fermions, which in our case is important 
for the top quarks in the ti production and SUSY particles in the cascades. To prop- 
erly account for the tau polarization, the generator was interfaced with the TAUOLA [29] 
package. 

We simulate again the four SUSY scenarios Scl-Sc4 described in section 2, with pro- 
duction cross sections times branching ratios as given in Fig. 5. For the SM background, 
again only the ti process was considered. Also the W + n-jet background will be impor- 
tant, but we have not included it since the standard Monte-Carlo generators do not give 
a realistic jet transverse energy spectrum for these processes, especially in the high trans- 
verse energy region. We are, however, convinced that any selection criterion that is able 
to suppress the ti background, will also successfully reject the W + n-jet background. 
Since the kinematical distributions of the SUSY cascades and the ti events are the same 
as in the neutral Higgs boson case, we can refer to Figs. 8-11 and the related discussion 
on the selection strategy to suppress the SM background. We will adopt here the same 
selection criteria. 

After eliminating the SM backgrounds, there still remains the more difficult task to 
discriminate the H ± bosons from the other particles in the SUSY cascade decays. Tau- 
leptons in the cascade backgrounds originate mostly from W and Z bosons. Taus coming 
from these particles will have a somewhat softer spectrum than the ones coming from 
a 170 GeV charged Higgs boson. Therefore we impose a lower limit on the transverse 
energy of the r-jet around the W mass, E^~ }et > 80 GeV. 

As a measure of the hardness of the charged hadron in the r-jet, we show in Fig. 16 the 
distribution of the fraction of the transverse energy E^~ }et that is carried by the one-prong 
charged particle in the jet. The full line corresponds to the SUSY signal [i.e. the cascade 
decays including the charged Higgs bosons] while the dashed line corresponds to the SUSY 
backgrounds [i.e. the same cascades where the decays of the heavier charginos/neutralinos 
to charged Higgs bosons and lighter charginos/neutralinos are forbidden]. Because of the 
polarization effects as explained above, this distribution peaks at high values for tau 
leptons coming from scalar charged Higgs bosons, while peaking at low values for tau 
leptons coming from vector-like W bosons. Requiring that at least 75% of the transverse 
energy of the r-jet be carried by the charged hadron, we can enhance the visibility of the 
signal with respect to the SUSY background. 

These considerations lead us to the following selection criteria in order to distinguish 
between signal and background events: 
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Figure 16: Distribution of the fraction of the i?J~ Jct that is carried by the charged particle 
in the jet. The full line represents the signal while the dashed line represents the SUSY 
background. The shaded areas contain the events left after applying the selection cuts, for 
an integrated luminosity of 30 fb~ l . The charged Higgs signal is clearly visible in Sc3 and 
Sc4, but not visible in Scl and Sc2. 



37 



i) The events should contain at least five jets. 
ii) The hardest jet in the event should have Et > 300 GeV. 

Hi) The transverse missing energy E™ ms in the event should be larger than 200 GeV. 

iv) We require exactly one hadronically decaying tau (1-prong), i.e. we demand a 
narrow jet within \q\ < 2.5, which should contain a hard charged track of pt > 5 
GeV in a cone of AR = 0.15 rad around the calorimeter jet axis, and it should be 
isolated, meaning no charged tracks with p T > 2 GeV are allowed in a cone of AR 
= 0.4 rad around the axis. 

v) The Et of the r-jet, defined as the Et reconstructed in a cone of AR = 0.4 rad 
around the jet axis, should be above 80 GeV. 

vi) More than 75% of the r-jet transverse energy should be carried by the charged 
track, i.e. we impose the cut p T /E T ~ >et > 0.75. 

Events that satisfy conditions can be efficiently triggered on using the jet and 

missing-energy triggers. 

In Fig. 16, the shaded areas represent the number of events left after all selection cuts, 
for the signal and the SUSY background. The charged Higgs boson signal is clearly visible 
in Sc3 and Sc4, but not in Scl and Sc2. This is because the production cross section for 
these scenarios are about 10 times smaller than for Sc3 and Sc4, as can be seen from 
Fig. 5. 

Clearly, the evidence for the charged Higgs boson is not as convincing as the invariant 
mass peaks we obtained in the neutral Higgs boson case. However, in the MSSM, the 
neutral and charged Higgs boson masses are not independent: since they are are connected 
through the relation M\± = M\ + [which is valid at leading order but is not strongly 
affected by radiative corrections], the observation of the neutral Higgs bosons can give 
definite information on the charged Higgs boson mass. If one then observes an excess in 
r-jet events after a selection as illustrated above, one can be rather confident that this 
excess originates from the production of charged Higgs particles 6 . 



6 Note that for relatively light charged Higgs bosons which can be produced through top quark decays, 
the situation can be slightly more complicated. Although the production rates of if ± bosons originating 
from cascade decays of squarks and gluinos can be enhanced as discussed in section 2.4, one has to deal 
with the additional events coming from direct production of tt pairs. More detailed analyses of this topic 
will be needed. 
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4. Conclusion 



We have analysed the detection prospects of the neutral and charged Higgs particles of 
the MSSM in single Higgs final states obtained through the cascade decays of the heavy 
squarks and gluinos that will be copiously produced at the LHC. We have discussed in 
detail the case of the production through the cascades involving the heavier charginos 
and neutralinos [originating from the decays of the strongly interacting SUSY particles] 
into the lighter charginos and neutralinos and a Higgs boson $. This can occur through 
the "big cascades", xl^xf ~ > Xi&xt + as well as through the "small cascades", 
X25 Xi Xi + 3") or through both cascades, depending on the considered scenario. For 
illustration, we have selected four scenarios which are expected to cover most of the 
situations which can occur in the MSSM: squarks either lighter or heavier than gluinos, 
and light charginos and neutralinos, either of the gaugino- or higgsino-type. 

We have shown that in these scenarios and with not too heavy Higgs particles, with 
masses M$ < 200-250 GeV, the cross sections times branching ratios for the production 
of at least one neutral or charged MSSM Higgs boson can be substantial at the LHC 7 , 
resulting in large numbers of signal events with an integrated luminosity of / L > 30 fb _1 , 
which is expected to be collected after a few years of running. 

We have performed a full Monte-Carlo simulation that takes into account the various 
signals as well as the SM and SUSY backgrounds and includes a fast simulation of some 
important aspects of the responce of the CMS detector at the LHC [similar conclusions 
are expected to be drawn for the ATLAS detector]. We have shown that these final 
states can indeed be detected in some representative MSSM scenarios. In particular, 
the heavier CP-even H and the pseudoscalar A bosons can be observed in the region of 
parameter space M4 ~ 150 GeV and tan/? ~ 5, which cannot be probed by searches 
using the standard Higgs production channels. Therefore, the search for Higgs bosons via 
the cascade decays could be complementary to the standard searches. 

This is exemplified in Fig. 17, in which we show the areas of the (tan /3, Ma) parameter 
space where the MSSM Higgs particles can be detected by the CMS collaboration with an 
integrated luminosity of J C ~ 100 fb _1 . Besides the usual areas where the Higgs bosons 
can be probed in SM-like production and decay processes, we exhibit the range of M4 
values where one can detect, in addition, the A, H and H ± bosons via the cascade decays 
in a specific MSSM scenario [namely Sc3, with the value of the parameter M 2 fixed to 
350 GeV]. The heavier neutral A, H particles can be probed for masses up to Ma ~ 220 

7 At the Tevatron, because of the reduced kinematical reach for squarks and gluinos, one would have 
access to these final states in much more limited areas of the parameter space. However, the production 
of the neutral Higgs boson through the little cascade, X2 ~ > Xi + ^7 might be kinematically possible and 
could be searched for with a high integrated luminosity. 
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Figure 17: The areas in the (M^,tan/?) parameter space where the MSSM Higgs bosons 
can be discovered at CMS with an integrated luminosity of 100 fb~ l . Various detection 
channels are shown in the case of the standard searches for the maximal stop mixing 
scenario. The right-hatched and cross-hatched regions show the areas where only the 
lightest h boson can be observed in these production channels. The left-hatched area is the 
region where the heavier CP-even H and pseudoscalar A bosons can be observed through 
the (big) cascade decays of squarks and gluinos in Sc3 with M2 = 350 GeV. A similar 
contour plot can be obtained in Sc4, where only the little cascade is at work. 
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GeV for the entire range of tan/5 values [as we have shown in our analysis, the cross 
sections times branching ratios depend only mildly on the value of this parameter], while 
the charged Higgs particle can be observed for masses up to M% ~ 200 GeV [the lighter 
h boson can of course be detected in the entire (tan (3, Ma) plane in this scenario]. As can 
be seen, the additional channels can fill some holes of the parameter space where these 
heavier Higgs particles cannot be observed through the standard processes. 

We have also discussed, though only at the level of production cross sections times 
decay branching fractions, MSSM Higgs boson production through the decays of heavier 
third-generation squark into the lighter ones, as well as the production of H ± bosons 
from the decays of top quarks originating from g,q cascade decays. While the rates for 
Higgs boson production from direct decays of third-generation squarks are in general 
smaller than in the previous case, they can be substantial for large tan/3 values in the 
case where top quarks originating from SUSY cascades decay into H ± bosons. These 
additional channels can therefore increase the yield for Higgs particles at the LHC We 
did not attempt to perform a Monte-Carlo simulation to verify whether the final states 
in these channels can indeed be detected in a realistic situation. 

The present study does not pretend to be exhaustive. To cover the entire possibilities 
for MSSM Higgs boson production through SUSY particle cascade decays, many more 
studies, in particular detailed experimental simulations, are required in order to fully 
cover this rather complicated subject. Here, we simply made a preliminary investigation 
in some selected scenarios, which indicates that the detection of Higgs particles in these 
processes is not at all hopeless and might even complement the searches through the 
intensively studied SM-like processes in some favourable cases. 

We stress again that these cascade processes are important not only because they 
represent a new source of Higgs bosons, but also because they will be very useful to 
measure the couplings of supersymmetric particles to Higgs bosons, which would be an 
essential ingredient to reconstruct the SUSY Lagrangian at low energies and to probe the 
theory at very high scales. 

In addition, if the branching ratios are sufficiently large, the decays of heavy neutrali- 
nos into bb pairs through Higgs boson intermediate states, xl,4 ~^ K H, A + xl — > bb + Xv 
contain information to reconstruct the masses of the heavier neutralinos and could as such 
help significantly in identifying the supersymmetric particle spectrum at the LHC. 

Hence, on their own merit, SUSY particle cascade processes involving Higgs bosons 
deserve detailed and dedicated studies in the future. 
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